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Abstract 

We construct finite mass, asymptotically flat black hole solutions in d = 5 Einstein- 
D \ Yang-Mills-Chern-Simons theory. Our results indicate the existence of a second order 

phase transition between Reissner-Nordstrom solutions and the non-Abelian black holes 
which generically are thermodynamically preferred. Some of the non-Abelian configu- 
rations are also stable under linear, spherically symmetric perturbations. In addition a 
■^j- ■ solution in closed form describing an extremal black hole with non-Abelian hair is found 

CN ■ for a special value of the Chern-Simons coupling constant. 
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Introduction.— The so called "no-hair" conjecture [1] states that an asymptotically flat, 
stationary black hole is uniquely described in terms of a small set of asymptotically measurable 
quantities. However, in recent years counterexamples to this conjecture were found in several 
theories, most of them containing non-Abelian matter fields. The first non-Abelian "hairy" 
black hole solutions within the framework of d — 4 SU(2) Einstein- Yang- Mills (EYM) theory, 
^ ■ were presented in [2]. Although these solutions were static and spherically symmetric with 
vanishing Yang-Mills (YM) charges, they were different from the Schwarzschild black hole and, 
therefore, not characterized exclusively by their total mass. However, all known asymptotically 
flat EYM solutions are perturbatively unstable and thus they do not contradict the spirit of 
the "no-hair" conjecture (see the review [3]). When considering instead a number d > 4 of 
spacetime dimensions, no finite energy asymptotically flat non-Abelian solutions are found 
[U |5] unless the action is supplemented with string-theory inspired higher order YM curvature 
terms [5J, in which case, again the solutions are classically unstable. 

In odd spacetime dimensions, the usual gauge field action can be augmented instead by a 
Chern-Simons (CS) term. Such (CS) terms appear in various supersymmetric theories. In the 
Abelian case this leads to some new features for rotating black holes only [TJ. In the non-Abelian 
case however, the CS term can affect the properties of solutions even in the static, spherically 
symmetric case. For the d = 5 case considered here, this allows the construction of finite 
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mass, asymptotically flat, non-Abelian black hole solutions which generically turn out to be 
thermodynamically favoured over the Reissner-Nordstrom Abelian configurations. Moreover, 
some of these solutions are stable against linear, spherically symmetric perturbations. 

The model.— In a 4 + 1 dimensional spacetime, the smallest simple gauge group supporting 
a nonvanishing CS term is SO (6). Then we consider a general EYMCS theory with Lagrangian 

C = —?—R*l - \*F AF- Kt h ... h (F hh A F hh A A hh (1) 
167rG 4 V 
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where A IJ are the £0(6) gauge fields, F IJ = dA IJ + gA IK A A KJ , G is gravitational constant, 
k the CS coefficient and g the gauge coupling constant. 

Our solutions are spherically symmetric with a line element 

ds 2 = -f dt 2 + hdr 2 + f 2 dnl (2) 

where /j are functions of r and t in general and dfl 2 = dip 2 + sin 2 ip(d9 2 + sin 2 9d<fi 2 ) is the line 
element of the three dimensional sphere. 

The general spherically symmetric Ansatz for the SO(6) YM field is given [8]. In the present 
work we restrict ourselves to a consistent truncation, 50(4) x 5*0(2), of that Ansatz. Apart 
from simplifying the picture, prominent new physical features are conveniently exposed in this 
truncation of the Ansatz, stated as, 

1 f w(r t) 4- 1 r? \ 
A = - v ' ' Eij—dx 1 + V(r, t)Y, 56 dt , with i, j = 1, . . . , 4 , (3) 
9 \ r r J 

Sy being the representation matrices of SO (4), and S 56 of the 5*0(2), subalgebras of 50(6). 
The Cartesian coordinates x % are related to the spherical coordinates (r, ip, 9, 0) as in flat space. 

Starting with static solutions, we have found it convenient to choose for the metric function 
in |2]), f x = 1/N(r), f 2 = r 2 , f = N(r)a 2 (r), with N(r) = 1 - m(r)/r 2 . Then the coupled 
static EYMCS equations of motion reduce to 



/ 1 i( (, T n (w 2 -l) 2 \ r-3 \ a' 3a V 2 



(raNw'Y = 2(TW{w2 1} + %kV'{w 2 - 1), (— )' = 2Ak{w 2 - l)w', 
r a 

where the prime indicates the derivative with respect to r and a 2 = 16nG/(3g 2 ) (a factor of 
l/g is also absorbed in k). The last equation above has the first integral 

V = ^(K + 8kw(w 2 -3)), (5) 

with K an integration constant. 

The set of equations (j3J) admit the scaling transformation 

r — V Xr, m — > X 2 m, a — >■ a, w — )■ w, V — > V/X, k — > Xk and a — > \a , (6) 
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which is used in what follows to set a = 1. 

We consider EYMCS black holes, with an event horizon located at r = > 0. In the 
vicinity of the event horizon, one finds the following expansion of the solutions 

m (r) =rl + m x (r - r h ) + . . . , a(r) =o h + -(r - r h ) + . . . , 

w(r) =w h + wi(r -r h ) + V(r) = v\{r - r h ) + . . . , (7) 

where v x = ^^Kzg)) ; mi = rfrW(i-^) a = ^si^)^ , while f or i arge va l U es 
of r, the expression of the solutions is 

m(r) = M - ^ + . . . , a(r) = 1 - ^ + • • • , «>(r) = ±1 + 4 + ■ ■ ■ > ^( r ) = -■? + •• • • 

In the above relations, at, Wh and J, M, Vq are parameters given by numerics which fix all 
higher order terms, while Q = K/2 =)= 8k. 

The only conserved quantities associated with these solutions are the mass M = %M and 
the electric charge Q = —Q, which is associated with the U(l) gauge symmetry generated by 
E 56 . Other quantities of interest are the chemical potential $ = — , the Hawking temperature 
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T H = ±a(r h )N'(r h ) and the entropy S = . 

The solutions.— The moduli space of black hole solutions includes the Reissner-Nordstrom 
black hole (hereafter RN), described by m(r) = M— V(r) = Vo—Q/r 2 , cr(r) = 1, w(r) = ±1. 

, i— v 1 /2 

This solution has an event horizon at rh — ( 4f + y ~ Q 2 ) > which becomes extremal for 

Q — We have found that for k > 1/8 and a given Q > 0, the RN black hole presents an 
instability with respect to static non-Abelian perturbations, for a critical value of the mass. 
This instability is found within the Ansatz (J3J) , for values of the magnetic gauge potential w(r) 
close to the value —1 everywhere, w(r) = — l + eW(r). The perturbation W(r) starts from some 
nonzero value at the horizon and vanishes at infinity, being a solution of the linear equation 

r(rNW')' - 4(1 - ^-)W = , (9) 

where N = 1 j-gk + The second term in ([9]) shows the existence of an effective mass 

term fi 2 for W near the horizon, with /i 2 ~ 1 — QnQ/r 2 ,. All physical solutions have fi 2 < 0, 
with Qjr\ = U(k) being a monotonic function of the CS coupling constant k. Then, for given 
K ^ 1/8, an instability occurs for a critical value of the mass to charge ratio of the RN solution, 
with M/Q = (1 + U 2 )/U. As k —> 1/8, one finds U = 1, while t/ ~ 1/4k for large «. No 
solutions of ([9]) are found for k < 1/8, or for perturbations of the form w(r) = +1 + eW(r), in 
which cases the effective mass for W is always real. 

This unstability signals the presence of a symmetry breaking branch of non-Abelian solutions 
bifurcating from the RN black hole. These nonperturbative configurations are found by solving 
numerically the Eqs. (j3j), using a shooting method. In this work we have restricted attention 
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Figure 1: The profiles of a typical non-Abelian (dotted line) solution with k = 0.5 is shown 
together with the corresponding Reissner- Nordstrom (RN) black hole (solid line) with the same 
mass and electric charge parameters M = 4.186 and Q = 2, respectively (the horizon radius 
of RN is smaller than that of the non-Abelian black hole). The intlet shows the corresponding 
potential for the perturbation equation (ITTj) below. 



to solutions with a monotonic behaviouiHJ of the magnetic gauge potential w(r). In contrast 
to other asymptotically flat non-Abelian black holes [2], some of the EYMCS solutions have 
no nodes in the magnetic gauge function w(r). A typical nodeless profile is shown in Fig. 1, 
together with the corresponding Abelian solution with the same mass and electric charge. 

Solutions smoothly interpolating between the asymptotics (jjTJ) and (jSJ) appear to exist for 
any value of the CS coupling constant k > 1/8. In a canonical ensemble, the non-Abelian black 
holes with a given k exist for a finite interval of {i.e. of the entropy) only. The detailed 
picture depends however on the ratio Q/k, with a critical value = 16k. For Q ^ Q^ c \ 
the temperature reaches its maximum at some intermediate value of the event horizon radius. 
Then a plot of the horizon function of the temperature reveals the existence of several 

branches of non-Abelian solutions. The upper branch ends in the critical RN solution with = 
^/Q/U(k). The lower branch possesses always a positive specific heat, the Hawking temperature 
vanishing there for a minimal value r[ mm ^ of the event horizon radius. As r — > r^™ m \ an 
extremal non-Abelian black hole solution with a regular horizon is approached. For Q < 
the near horizon expansion of the solutions implies r[ mm ' ) = 4a/3k(1 — w\) / ^/QAk 2 — w 2 hl where 
w h satisfies the equation (64k 2 — Q 2 + 2 k(1 + w h ) 2 (128 k 2 (w h — 2)+ w^w 2 ,— 2w h + 3)) = 0; for 
Q > Q^ c \ the limiting solution has rj l " Mn ' ) = y/Q — 16k and Wh = 1- Some of these features are 
shown in Fig. 2 where we plot the reduced area of the horizon = 2n 2 rf l /Q 3 ^ 2 as a function 
of the dimensionless temperature t H for a fixed CS coupling constant and several values of 

1 Although there exist solutions where w has local extrema, it is likely that these are always thermodynami- 
cally disfavoured because spatial oscillations in w increase the total mass. 
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Figure 2: The dimensionless horizon area is plotted vs. the scaled temperature for the non- 
Abelian solutions with several values of the electric charge and a given value of k. The branch 
of Reissner- Nordstrom solutions is also shown. 



the charge parameter Q (both an and tn are invariant under the scaling symmetry ((6])). The 
branch of RN solutions is also shown there. 

Furthermore, it turns out that the free energy F = M. — TnS of a RN solution is larger 
than the free energy of a lower branch non-Abelian solution with the same temperature and 
electric charge, except for configurations with k close to 1/8 and small enough values of the 
charge, Q < /3. Therefore the non-Abelian black holes are generically preferred. This 
reveals the existence of a second order phase transition between the RN solutions and the non- 
Abelian solutions. These aspects are exhibited in Fig. 3 where the dimensionless free energy 
/ = GF/Q is plotted as a function of the dimensionless temperature t# = T H y/Q for two values 
of k. Moreover, as seen e.g. in Fig. 1, for the same values of the mass and electric charge, 
the RN solution typically has a smaller event horizon radius (and thus a smaller entropy), 
than the non-Abelian black hole. The parameter J which enters the large r asymptotics of the 
magnetic gauge potential w(r) increases from zero (for the critical RN solution) to a maximal 
value approached at T# = 0. 

The overall picture is somehow different for Q = Q^ c \ in which case, despite the presence 
of an electric charge, the non-Abelian black holes behave in a similar way to the vacuum 
Schwarzschild-Tangherlini solution, with a single branch of thermally unstable configurations, 
see Fig. 2. In the limit — > 0, these black holes approach a set of globally regular particle-like 
solutions, whose mass is an almost linear function of k. 

In the special case k = 1/8, following the approach in p], one finds the following exact 

2 Note that the non-Abelian solutions with large enough temperature have no RN counterparts, see Fig. 3. 
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Figure 3: The scaled free energy is plotted vs. the scaled temperature for the Reissner- 
Nordstrom solutions (dotted line) and the two non-Abelian sets of solutions with different 
values of the Chern-Simons coupling constant k and the same charge parameter Q = 3.8. 



solution of the EYMCS equations within the general Ansatz (J2J), ©: 

where Q > 2 and J are arbitrary parameters. This describes an extremal black hole with 
non-Abelian hair, the regular event horizon being at r = (in these coordinates). The mass, 
electric charge and entropy of this non-Abelian deformation of the extremal RN solution are 
M = Q/g = ^tQ, S = 8it 3 (Q — 2) 3 / 2 /3g 2 (note that J does not enter any physical quantity). 
In the limit Q — > 2, the solution (|10|) describes a particle-like soliton with a regular origin. 

The stability of solutions.— The fact that we have found nodeless solutions suggest 
the existence of non-Abelian configurations stable against spherically symmetric perturbations. 
In examining such time-dependent fluctuations, we consider the metric Ansatz (|2j) with j\ = 
1 - '^A h = r 2 , fo = (1 - !Z $r i )a 2 (r,t), the YM Ansatz © and the following perturbed 
variables^ 



m(r,t) = m(r) + em 1 (r)e int , a(r,t) = <r(r)(l + ea 1 (r)e lUt ) 



w(r,t) = w(r) + ew 1 (r)e l ' l \ V(r,t) = V(r) + eV^e 1 ^. 

One finds that in the system of linearised EYMCS equations, the functions m 1; V\ and o\ can 
be eliminated in favor of wi(r), leading to a single Schrodinger equation for w\\ 

-^ + U n (p)x = n 2 X , (11) 



3 The corresponding problem for the SU(2) EYM system has been considered in [5]. 
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where x = w w/r, dr /dp = Na and 




( 9 -w' 2 + 192k 2 - |) 



r 2 V' 2 (l-6w' 2 )- 
A& 2 . 



(12) 



a 



The potential above is regular in the entire range — oo < p < oo. Near the event horizon, one 
finds Uq — > 0; for large values of p the potential is positive and bounded. Standard results from 
quantum mechanics [TD] further imply that there are no negative eigenvalues for Q 2 (and then 
no unstable modes) if the potential Uq is everywhere positive. 

Although the potential (f!2|) is not positive definite for all values of Q, k, we have found 
numerically that the condition Un > is fulfilled by some of the nodeless solutions (see the Fig. 
1 for a such a configuration). Therefore at least some of our solutions are linearly stable. The 
full picture is, however, quite complicated and a detailed discussion will be presented elsewhere. 
We note only that, by using the approach in [TT] we have found that the EYMCS solutions 
with one node in w(r) we have constructed are indeed unstable. 

Further remarks.— The main purpose of this work was to provide an example of stable, 
non-Abelian black holes without scalar fields, in asymptotically Minkowski spacetime. 

One should note that, despite the different asymptotic structure of spacetime and the dif- 
ferent horizon topology, the solutions in this work have some similarities with the colorful black 
holes with charge in Anti-de Sitter (AdS) space [T2"| [T5] . which provide a model of holographic 
superconductors. In both cases, an Abelian gauge symmetry is spontaneously broken near 
a black hole horizon with the appearance of a condensate of non-Abelian gauge fields there, 
leading to a second order phase transition. It remains an interesting open problem to clarify 
if the asymptotically flat EYMCS black holes may also provide useful analogies to phenomena 
observed in condensed matter physics. However, given the presence of several branches, the 
picture is more complicated for asymptotically flat solutions and we could not find so far simple 
universal relations between the relevant parameters as those found in [12] in the AdS case. 

The black holes in this work admit a straightforward generalisation with a cosmological 
constant. For A < 0, the basic properties of asymptotically AdS configurations were discussed 
in [S] . In this case the EYMCS model with a special value of k can be thought of as a truncation 
of the Af = 8, d = 5 gauged supergravity [HI US], with all scalars there taking constant values. 
We have found that most of the properties of the asymptotically flat solutions, in particular the 
existence of stable configurations, hold also in the AdS case. In this case the solutions within 
the restricted SO (4) x 0(2) Ansatz ([3]) contain already all relevant features of the full SO (6) 
configurations. 

While the existence of hairy non-Abelian black holes with AdS asymptotics is not a surprise 
in view of the results in [5], [16], we have found that, the A = solutions presented here admit 
generalisations with de-Sitter asymptotics (A > 0) as well. These black holes, possesing a 
regular cosmological event horizon at r = r c > r^, were constructed within the same Ansatz as 
in the asymptotically flat case. By using an inflationary coordinate system, one finds also the 
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following generalisation of fllUp . with k = 1/8 again and 



d S 2 = -^—Adr 2 + r 2 dn 2 3 ) - F 2 (r, t)dt\ w(r) = J — — , V(r, t) = F(r, t), (13) 

where F(r, t) _1 = 1 + e~ 2m + pqr^p ) and A = 6// 2 . This configuration describes a 

non-Abelian deformation of the extremal RN-de Sitter black holes in [T7j. By using the methods 
in [T7] , [18] , one finds that ( TT3l) shares the basic features of the J = Abelian configuration. 

Finally, we conjecture that similar EYMCS solutions exist in higher (odd) dimensions, the 
features encountered here being universal. 
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